svyset _n [pweight=weight], vce(linearized) singleunit(missing)
gen Media_Bias1=TAM152
gen Media_Bias2=TAM215
gen bias_change=Media_Bias2 - Media_Bias1
recode pid7 (8=.) (98=.) (99=.), gen(PID)
recode newsint (4=0) (3=1) (2=2) (1=3), gen(polint)
gen pidXpolint=PID*polint

gen sdem=(PID==1)
gen wdem=(PID==2)
gen ldem=(PID==3)
gen ind=(PID==4)
gen lrep=(PID==5)
gen wrep=(PID==6)
gen srep=(PID==7)

**Models**

svy linearized, subpop(if tookpost==1)  : tobit Media_Bias1 sdem wdem ldem lrep wrep srep polint daysprior
outreg using mbiasa1, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) replace
svy linearized, subpop(if tookpost==1)  : tobit Media_Bias2 sdem wdem ldem lrep wrep srep polint dayspost
outreg using mbiasa1, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1 & polint==3)  : tobit Media_Bias1 sdem wdem ldem lrep wrep srep daysprior
outreg using mbiasa1, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1 & polint==3)  : tobit Media_Bias2 sdem wdem ldem lrep wrep srep dayspost
outreg using mbiasa1, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1 & polint<3)  : tobit Media_Bias1 sdem wdem ldem lrep wrep srep daysprior
outreg using mbiasa1, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge 
svy linearized, subpop(if tookpost==1 & polint<3)  : tobit Media_Bias2 sdem wdem ldem lrep wrep srep dayspost
outreg using mbiasa1, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge


svy linearized, subpop(if tookpost==1)  : tobit Media_Bias2 Media_Bias1 sdem wdem ldem lrep wrep srep polint daysprior dayspost
outreg using mbiasc1, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) replace
svy linearized, subpop(if tookpost==1 & polint==3)  : tobit Media_Bias2 Media_Bias1 sdem wdem ldem lrep wrep srep daysprior dayspost
outreg using mbiasc1, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1 & polint<3)  : tobit Media_Bias2 Media_Bias1 sdem wdem ldem lrep wrep srep daysprior dayspost
outreg using mbiasc1, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge

svy linearized, subpop(if tookpost==1)  : regress bias_change sdem wdem ldem lrep wrep srep days
outreg using mbiasb1, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) replace
svy linearized, subpop(if tookpost==1 & polint==3)  : regress bias_change sdem wdem ldem lrep wrep srep days
outreg using mbiasb1, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1 & polint<3)  : regress bias_change sdem wdem ldem lrep wrep srep days
outreg using mbiasb1, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge



*Models with Covariates**

svy linearized, subpop(if tookpost==1)  : regress Media_Bias1 sdem wdem ldem lrep wrep srep polint age educ gender2 race_black hispanic2
outreg using mbias1, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) replace
svy linearized, subpop(if tookpost==1)  : regress Media_Bias2 sdem wdem ldem lrep wrep srep polint age educ gender2 race_black hispanic2
outreg using mbias1, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1 & polint==3)  : regress Media_Bias1 sdem wdem ldem lrep wrep srep age educ gender2 race_black hispanic2
outreg using mbias1, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1 & polint==3)  : regress Media_Bias2 sdem wdem ldem lrep wrep srep age educ gender2 race_black hispanic2
outreg using mbias1, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1 & polint<3)  : regress Media_Bias1 sdem wdem ldem lrep wrep srep age educ gender2 race_black hispanic2
outreg using mbias1, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1 & polint<3)  : regress Media_Bias2 sdem wdem ldem lrep wrep srep age educ gender2 race_black hispanic2
outreg using mbias1, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge


svy linearized, subpop(if tookpost==1)  : regress Media_Bias2 Media_Bias1 sdem wdem ldem lrep wrep srep polint age educ gender2 race_black hispanic2
outreg using mbias2, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) replace
svy linearized, subpop(if tookpost==1)  : regress bias_change sdem wdem ldem lrep wrep srep polint age educ gender2 race_black hispanic2
outreg using mbias2, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1  & polint==3)  : regress Media_Bias2 Media_Bias1 sdem wdem ldem lrep wrep srep polint age educ gender2 race_black hispanic2
outreg using mbias2, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1  & polint==3)  : regress bias_change sdem wdem ldem lrep wrep srep polint age educ gender2 race_black hispanic2
outreg using mbias2, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1 & polint<3)  : regress Media_Bias2 Media_Bias1 sdem wdem ldem lrep wrep srep age educ gender2 race_black hispanic2
outreg using mbias2, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1 & polint<3)  : regress bias_change sdem wdem ldem lrep wrep srep age educ gender2 race_black hispanic2
outreg using mbias2, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge



*testing interaction effect for PID and polint**

svy linearized, subpop(if tookpost==1)  : regress Media_Bias1 pid7 polint pidXpolint  age educ gender2 race_black hispanic2
outreg using mbias3, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) replace
svy linearized, subpop(if tookpost==1)  : regress Media_Bias2 pid7 polint pidXpolint  age educ gender2 race_black hispanic2
outreg using mbias3, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1)  : regress Media_Bias2 Media_Bias1 pid7 polint pidXpolint  age educ gender2 race_black hispanic2
outreg using mbias3, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1)  : regress bias_change pid7 polint pidXpolint  age educ gender2 race_black hispanic2
outreg using mbias3, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge

svy linearized, subpop(if tookpost==1)  : regress Media_Bias1 sdem wdem ldem lrep wrep srep polint intXsdem intXwdem intXldem intXlrep intXwrep intXsrep
outreg using mbias4, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) replace
svy linearized, subpop(if tookpost==1)  : regress Media_Bias2 sdem wdem ldem lrep wrep srep polint intXsdem intXwdem intXldem intXlrep intXwrep intXsrep
outreg using mbias4, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1)  : regress Media_Bias2 Media_Bias1 sdem wdem ldem lrep wrep srep polint intXsdem intXwdem intXldem intXlrep intXwrep intXsrep
outreg using mbias4, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1)  : regress bias_change sdem wdem ldem lrep wrep srep polint intXsdem intXwdem intXldem intXlrep intXwrep intXsrep
outreg using mbias4, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge

**Alternative specification with two category interest**
gen int2Xsdem=polint2*sdem
gen int2Xwdem=polint2*wdem
gen int2Xldem=polint2*ldem
gen int2Xlrep=polint2*lrep
gen int2Xwrep=polint2*wrep
gen int2Xsrep=polint2*srep

svy linearized, subpop(if tookpost==1)  : regress Media_Bias1 sdem wdem ldem lrep wrep srep polint int2Xsdem int2Xwdem int2Xldem int2Xlrep int2Xwrep int2Xsrep
outreg using mbias5, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) replace
svy linearized, subpop(if tookpost==1)  : regress Media_Bias2 sdem wdem ldem lrep wrep srep polint int2Xsdem int2Xwdem int2Xldem int2Xlrep int2Xwrep int2Xsrep
outreg using mbias5, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1)  : regress Media_Bias2 Media_Bias1 sdem wdem ldem lrep wrep srep polint int2Xsdem int2Xwdem int2Xldem int2Xlrep int2Xwrep int2Xsrep
outreg using mbias5, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1)  : regress bias_change sdem wdem ldem lrep wrep srep polint int2Xsdem int2Xwdem int2Xldem int2Xlrep int2Xwrep int2Xsrep
outreg using mbias5, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge

**models accounting for time**

svy linearized, subpop(if tookpost==1)  : regress Media_Bias1 sdem wdem ldem lrep wrep srep polint daysprior 
svy linearized, subpop(if tookpost==1)  : regress Media_Bias2 sdem wdem ldem lrep wrep srep polint dayspost

svy linearized, subpop(if tookpost==1 & polint==3)  : regress Media_Bias1 sdem wdem ldem lrep wrep srep daysprior 
svy linearized, subpop(if tookpost==1 & polint==3)  : regress Media_Bias2 sdem wdem ldem lrep wrep srep dayspost

svy linearized, subpop(if tookpost==1 & polint<3)  : regress Media_Bias1 sdem wdem ldem lrep wrep srep daysprior  
svy linearized, subpop(if tookpost==1 & polint<3)  : regress Media_Bias2 sdem wdem ldem lrep wrep srep dayspost

svy linearized, subpop(if tookpost==1)  : regress Media_Bias2 Media_Bias1 sdem wdem ldem lrep wrep srep polint days
svy linearized, subpop(if tookpost==1 & polint==3)  : regress Media_Bias2 Media_Bias1 sdem wdem ldem lrep wrep srep days
svy linearized, subpop(if tookpost==1 & polint<3)  : regress Media_Bias2 Media_Bias1 sdem wdem ldem lrep wrep srep days

svy linearized, subpop(if tookpost==1)  : regress bias_change sdem wdem ldem lrep wrep srep polint days
svy linearized, subpop(if tookpost==1 & polint==3)  : regress bias_change sdem wdem ldem lrep wrep srep days
svy linearized, subpop(if tookpost==1 & polint<3)  : regress bias_change sdem wdem ldem lrep wrep srep daysprior days

***Tobit Models***
svy linearized, subpop(if tookpost==1)  : tobit Media_Bias1 sdem wdem ldem lrep wrep srep polint, ul ll 
outreg using mbiasT1, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) replace 
svy linearized, subpop(if tookpost==1)  : tobit Media_Bias2 sdem wdem ldem lrep wrep srep polint, ul ll
outreg using mbiasT1, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge

svy linearized, subpop(if tookpost==1 & polint==3)  : tobit Media_Bias1 sdem wdem ldem lrep wrep srep, ul ll
outreg using mbiasT1, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge  
svy linearized, subpop(if tookpost==1 & polint==3)  : tobit Media_Bias2 sdem wdem ldem lrep wrep srep, ul ll
outreg using mbiasT1, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge

svy linearized, subpop(if tookpost==1 & polint<3)  : tobit Media_Bias1 sdem wdem ldem lrep wrep srep, ul ll
outreg using mbiasT1, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge   
svy linearized, subpop(if tookpost==1 & polint<3)  : tobit Media_Bias2 sdem wdem ldem lrep wrep srep, ul ll 
outreg using mbiasT1, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge

svy linearized, subpop(if tookpost==1)  : tobit Media_Bias2 Media_Bias1 sdem wdem ldem lrep wrep srep polint, ul ll 
outreg using mbiasT2, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) replace 
svy linearized, subpop(if tookpost==1 & polint==3)  : tobit Media_Bias2 Media_Bias1 sdem wdem ldem lrep wrep srep, ul ll 
outreg using mbiasT2, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1 & polint<3)  : tobit Media_Bias2 Media_Bias1 sdem wdem ldem lrep wrep srep, ul ll 
outreg using mbiasT2, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge

svy linearized, subpop(if tookpost==1)  : tobit bias_change sdem wdem ldem lrep wrep srep polint, ul ll
outreg using mbiasT2, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1 & polint==3)  : tobit bias_change sdem wdem ldem lrep wrep srep, ul ll 
outreg using mbiasT2, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1 & polint<3)  : tobit bias_change sdem wdem ldem lrep wrep srep, ul ll  
outreg using mbiasT2, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge

***Tobit Models with coefplot***

svy linearized, subpop(if tookpost==1)  : tobit Media_Bias1 sdem wdem ldem lrep wrep srep polint, ul ll 
estimates store BiasT1A
svy linearized, subpop(if tookpost==1)  : tobit Media_Bias2 sdem wdem ldem lrep wrep srep polint, ul ll
estimates store BiasT2A

svy linearized, subpop(if tookpost==1 & polint==3)  : tobit Media_Bias1 sdem wdem ldem lrep wrep srep, ul ll
estimates store BiasT1B 
svy linearized, subpop(if tookpost==1 & polint==3)  : tobit Media_Bias2 sdem wdem ldem lrep wrep srep, ul ll
estimates store BiasT2B

svy linearized, subpop(if tookpost==1 & polint<3)  : tobit Media_Bias1 sdem wdem ldem lrep wrep srep, ul ll
estimates store BiasT1C
svy linearized, subpop(if tookpost==1 & polint<3)  : tobit Media_Bias2 sdem wdem ldem lrep wrep srep, ul ll 
estimates store BiasT2C
coefplot (BiasT1A, label(Bias t1)) (BiasT2A, label(Bias t2)), bylabel(All Respondents) drop(_cons) xline(0) || BiasT1B BiasT2B, bylabel(High Interest) drop(_cons) xline(0)|| BiasT1C BiasT2C, bylabel(Low Interest) drop(_cons) xline(0)

coefplot (BiasT1A, label(Bias t1)) (BiasT2A, label(Bias t2)), bylabel(All Respondents) drop(_cons) xline(0) 
coefplot (BiasT1B, label(Bias t1)) (BiasT2B, label(Bias t2)), bylabel(All Respondents) drop(_cons) xline(0) 
coefplot (BiasT1C, label(Bias t1)) (BiasT2C, label(Bias T2)), bylabel(High Interest) drop(_cons) xline(0)

***Final coefplots**
coefplot (BiasT1A, label(Bias t1)), bylabel(Media Bias t1) || (BiasT2A, label(Bias t2)), bylabel(Media Bias t2) drop(_cons) xline(0) graphregion(color(white))
coefplot (BiasT1B, label(Bias t1)), bylabel(Media Bias t1) || (BiasT2B, label(Bias t2)) || (BiasT1C, label (Bias t1)) || (BiasT2C, label (Bias t2)) bylabel(Media Bias t2) drop(_cons) xline(0) graphregion(color(white))
coefplot (BiasT1C, label(Bias t1)), bylabel(Media Bias t1) || (BiasT2C, label(Bias t2)), bylabel(Media Bias t2) drop(_cons) xline(0) graphregion(color(white))



**Change Models Coefplot**
svy linearized, subpop(if tookpost==1)  : tobit Media_Bias2 Media_Bias1 sdem wdem ldem lrep wrep srep polint, ul ll 
estimates store L1
svy linearized, subpop(if tookpost==1 & polint==3)  : tobit Media_Bias2 Media_Bias1 sdem wdem ldem lrep wrep srep, ul ll 
estimates store L2
svy linearized, subpop(if tookpost==1 & polint<3)  : tobit Media_Bias2 Media_Bias1 sdem wdem ldem lrep wrep srep, ul ll 
estimates store L3

coefplot L2, label(High Interest)), bylabel(High Interest) || (L3, label(Low Interest)), bylabel(Low Interest) drop(_cons) xline(0) graphregion(color(white))

svy linearized, subpop(if tookpost==1)  : tobit bias_change sdem wdem ldem lrep wrep srep polint, ul ll
estimates store C1B
svy linearized, subpop(if tookpost==1 & polint==3)  : tobit bias_change sdem wdem ldem lrep wrep srep, ul ll 
estimates store C2B
svy linearized, subpop(if tookpost==1 & polint<3)  : tobit bias_change sdem wdem ldem lrep wrep srep, ul ll  
estimates store C3B

coefplot (C2B, label(High Interest)), bylabel(High Interest) || (C3B, label(Low Interest)), bylabel(Low Interest) drop(_cons) xline(0) graphregion(color(white)) 

coefplot L2, bylabel(Highly Interested) || L3, bylabel(Highly Interested) drop(_cons) xline(0) graphregion(color(white)) 
coefplot C2B, bylabel(Less Interested) || C3B, bylabel(Less Interested) drop(_cons) xline(0) graphregion(color(white)) 

**Tobit Interactions***






svy linearized, subpop(if tookpost==1)  : tobit Media_Bias1 sdem wdem ldem lrep wrep srep polint intXsdem intXwdem intXldem intXlrep intXwrep intXsrep, ul ll 
outreg using mbiasT3, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) replace
svy linearized, subpop(if tookpost==1)  : tobit Media_Bias2 sdem wdem ldem lrep wrep srep polint intXsdem intXwdem intXldem intXlrep intXwrep intXsrep, ul ll 
outreg using mbiasT3, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1)  : tobit Media_Bias2 Media_Bias1 sdem wdem ldem lrep wrep srep polint intXsdem intXwdem intXldem intXlrep intXwrep intXsrep, ul ll 
outreg using mbiasT3, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1)  : tobit bias_change sdem wdem ldem lrep wrep srep polint intXsdem intXwdem intXldem intXlrep intXwrep intXsrep, ul ll 
outreg using mbiasT3, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge


**Tobit test for days**

svy linearized, subpop(if tookpost==1)  : tobit Media_Bias1 sdem wdem ldem lrep wrep srep polint daysprior, ul ll 
outreg using mbiasT1, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) replace 
svy linearized, subpop(if tookpost==1)  : tobit Media_Bias2 sdem wdem ldem lrep wrep srep polint days, ul ll
outreg using mbiasT1, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge

svy linearized, subpop(if tookpost==1 & polint==3)  : tobit Media_Bias1 sdem wdem ldem lrep wrep srep daysprior, ul ll
outreg using mbiasT1, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge  
svy linearized, subpop(if tookpost==1 & polint==3)  : tobit Media_Bias2 sdem wdem ldem lrep wrep srep days, ul ll
outreg using mbiasT1, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge

svy linearized, subpop(if tookpost==1 & polint<3)  : tobit Media_Bias1 sdem wdem ldem lrep wrep srep daysprior, ul ll
outreg using mbiasT1, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge   
svy linearized, subpop(if tookpost==1 & polint<3)  : tobit Media_Bias2 sdem wdem ldem lrep wrep srep days, ul ll 
outreg using mbiasT1, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge


svy linearized, subpop(if tookpost==1)  : tobit Media_Bias2 Media_Bias1 sdem wdem ldem lrep wrep srep polint days, ul ll 
outreg using mbiasT2, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) replace 
svy linearized, subpop(if tookpost==1 & polint==3)  : tobit Media_Bias2 Media_Bias1 sdem wdem ldem lrep wrep srep days, ul ll 
outreg using mbiasT2, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1 & polint<3)  : tobit Media_Bias2 Media_Bias1 sdem wdem ldem lrep wrep srep days, ul ll 
outreg using mbiasT2, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge

svy linearized, subpop(if tookpost==1)  : tobit bias_change sdem wdem ldem lrep wrep srep polint days, ul ll
outreg using mbiasT2, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1 & polint==3)  : tobit bias_change sdem wdem ldem lrep wrep srep days, ul ll 
outreg using mbiasT2, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1 & polint<3)  : tobit bias_change sdem wdem ldem lrep wrep srep days, ul ll  
outreg using mbiasT2, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge


svy linearized, subpop(if tookpost==1)  : regress bias_change sdem wdem ldem lrep wrep srep polint
estimates store C1
svy linearized, subpop(if tookpost==1 & polint==3)  : regress bias_change sdem wdem ldem lrep wrep srep
estimates store C2
svy linearized, subpop(if tookpost==1 & polint<3)  : regress bias_change sdem wdem ldem lrep wrep srep
estimates store C3

coefplot (C2, label(High Interest)), bylabel(High Interest) || (C3, label(Low Interest)), bylabel(Low Interest) drop(_cons) xline(0) 


*Tobit Models with Covariates**

svy linearized, subpop(if tookpost==1)  : tobit Media_Bias1 sdem wdem ldem lrep wrep srep polint age educ gender2 race_black hispanic2, ul ll
outreg using mbias1A, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) replace
svy linearized, subpop(if tookpost==1)  : tobit Media_Bias2 sdem wdem ldem lrep wrep srep polint age educ gender2 race_black hispanic2, ul ll
outreg using mbias1A, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1 & polint==3)  : tobit Media_Bias1 sdem wdem ldem lrep wrep srep age educ gender2 race_black hispanic2, ul ll
outreg using mbias1A, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1 & polint==3)  : tobit Media_Bias2 sdem wdem ldem lrep wrep srep age educ gender2 race_black hispanic2, ul ll
outreg using mbias1A, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1 & polint<3)  : tobit Media_Bias1 sdem wdem ldem lrep wrep srep age educ gender2 race_black hispanic2, ul ll
outreg using mbias1A, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1 & polint<3)  : tobit Media_Bias2 sdem wdem ldem lrep wrep srep age educ gender2 race_black hispanic2, ul ll
outreg using mbias1A, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge


svy linearized, subpop(if tookpost==1)  : tobit Media_Bias2 Media_Bias1 sdem wdem ldem lrep wrep srep polint age educ gender2 race_black hispanic2, ul ll
outreg using mbias2, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) replace
svy linearized, subpop(if tookpost==1)  : tobit bias_change sdem wdem ldem lrep wrep srep polint age educ gender2 race_black hispanic2, ul ll
outreg using mbias2, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1  & polint==3)  : tobit Media_Bias2 Media_Bias1 sdem wdem ldem lrep wrep srep polint age educ gender2 race_black hispanic2, ul ll
outreg using mbias2, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1  & polint==3)  : tobit bias_change sdem wdem ldem lrep wrep srep polint age educ gender2 race_black hispanic2, ul ll
outreg using mbias2, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1 & polint<3)  : tobit Media_Bias2 Media_Bias1 sdem wdem ldem lrep wrep srep age educ gender2 race_black hispanic2, ul ll
outreg using mbias2, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1 & polint<3)  : tobit bias_change sdem wdem ldem lrep wrep srep age educ gender2 race_black hispanic2, ul ll
outreg using mbias2, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge

svy linearized, subpop(if tookpost==1)  : regress bias_change sdem wdem ldem lrep wrep srep polint age educ gender2 race_black hispanic2
outreg using mbias9, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) replace
svy linearized, subpop(if tookpost==1  & polint==3)  : regress bias_change sdem wdem ldem lrep wrep srep polint age educ gender2 race_black hispanic2
outreg using mbias9, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1 & polint<3)  : regress bias_change sdem wdem ldem lrep wrep srep age educ gender2 race_black hispanic2
outreg using mbias9, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge


svy linearized, subpop(if tookpost==1)  : tobit Media_Bias1 pid7c polint pidXpolint  age educ gender2 race_black hispanic2, ul ll
outreg using mbias53, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) replace
svy linearized, subpop(if tookpost==1)  : tobit Media_Bias2 pid7c polint pidXpolint  age educ gender2 race_black hispanic2, ul ll
outreg using mbias53, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1)  : tobit Media_Bias2 Media_Bias1 pid7c polint pidXpolint  age educ gender2 race_black hispanic2, ul ll
outreg using mbias53, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1)  : regress bias_change pid7c polint pidXpolint  age educ gender2 race_black hispanic2
outreg using mbias63, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) replace


***Trump Vote Models***
svy linearized, subpop(if tookpost==1)  : tobit Media_Bias1 HillaryVote polint age educ gender2 race_black hispanic2, ul ll
outreg using mbias31, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) replace
svy linearized, subpop(if tookpost==1)  : tobit Media_Bias2 HillaryVote polint age educ gender2 race_black hispanic2, ul ll
outreg using mbias31, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1 & polint==3)  : tobit Media_Bias1 HillaryVote age educ gender2 race_black hispanic2, ul ll
outreg using mbias31, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1 & polint==3)  : tobit Media_Bias2 HillaryVote age educ gender2 race_black hispanic2, ul ll
outreg using mbias31, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1 & polint<3)  : tobit Media_Bias1 HillaryVote age educ gender2 race_black hispanic2, ul ll
outreg using mbias31, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1 & polint<3)  : tobit Media_Bias2 HillaryVote age educ gender2 race_black hispanic2, ul ll
outreg using mbias31, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge


svy linearized, subpop(if tookpost==1)  : tobit Media_Bias2 Media_Bias1 HillaryVote polint age educ gender2 race_black hispanic2
outreg using mbias32, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) replace
svy linearized, subpop(if tookpost==1 & polint==3)  : tobit Media_Bias2 Media_Bias1 HillaryVote age educ gender2 race_black hispanic2
outreg using mbias32, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1 & polint<3)  : tobit Media_Bias2 Media_Bias1 HillaryVote age educ gender2 race_black hispanic2
outreg using mbias32, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge


svy linearized, subpop(if tookpost==1)  : regress bias_change HillaryVote polint age educ gender2 race_black hispanic2
outreg using mbias33, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) replace
svy linearized, subpop(if tookpost==1 & polint==3)  : regress bias_change HillaryVote age educ gender2 race_black hispanic2
outreg using mbias33, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1 & polint<3)  : regress bias_change HillaryVote age educ gender2 race_black hispanic2
outreg using mbias33, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge


svy linearized, subpop(if tookpost==1)  : tobit Media_Bias1 HillaryVote polint IntXHillary age educ gender2 race_black hispanic2
outreg using mbias35, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) replace
svy linearized, subpop(if tookpost==1)  : tobit Media_Bias2 HillaryVote polint  IntXHillary age educ gender2 race_black hispanic2
outreg using mbias35, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1)  : tobit Media_Bias2 Media_Bias1 HillaryVote polint IntXHillary age educ gender2 race_black hispanic2
outreg using mbias35, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) replace
svy linearized, subpop(if tookpost==1)  : regress bias_change HillaryVote polint IntXHillary age educ gender2 race_black hispanic2
outreg using mbias35, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) replace




svy linearized, subpop(if tookpost==1)  : tobit Media_Bias1 HillaryVote polint IntXHillary age educ gender2 race_black hispanic2
outreg using mbias31, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge
svy linearized, subpop(if tookpost==1)  : tobit Media_Bias2 HillaryVote polint IntXHillary age educ gender2 race_black hispanic2
outreg using mbias31, se bdec(2 2 2 2) summstat(r2) summtitle(R-squared) summdec(2 2) merge